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This work is devoted to the characterization of the SI:H system obtained by high-fluence, low-energy, 
hydrogen implantation Into single-crystal silicon. The implanted hydrogen profile and the ones resulting 
after thermal annealing in the range 100-300 'C are detected by secondary-ion mass spectrometry and 
elastic-recoil detection analysis. The displacement field in the crystal, measured by channeling 
Rutherford-backscattering spectrometry, is found to depend on the direct radiation damage, the extend- 
ed defects formed alter ion implantation (revealed by transmission electron microscopy), and the im- 
planted species. The contribution to the displacement field due to hydrogen-related defects has a charac- 
teristic "reverse annealing" in the range 100-400 *C, essentially due to their formation kinetics. 



t INTRODUCTION 

Though "none of the numerous hydrogen-related vi- 
bration bonds have yet been assigned to a specific de- 
fect", 1 there is now a reasonable consensus on hydrogen 
configurations in single-crystal silicon tX low concentra- 
tion. 

(1) The bond-centered CBC) configuration, typical of in- 
trinsic or /Mype silicon [Fig. 1(a) and Refs. 2 and 3]. 

(2) The antibonding (AB) configuration, where hydro- 
gen is kept close to an interstitial tetrahedral {T d ) posi- 
tion by an antibonding state of silicon [Fig. 1(b) and Ref. 
4]. 

(3) The HJ configuration, resulting from pairing of 
close, ionized, AB (H~) and BC (H + ) hydrogen atoms 
[Fig. 1(c) and Ref. 3]. 

(4) The aflame Si-H configuration, resulting from hy- 
drogen reaction with a silicon dangling bond [Fig. 1(d)]. 

The standard theory of shallow dopants ascribes accep- 
tor passivation to a positively ionized BC hydrogen in the 
vicinity of a negatively ionized acceptor [Fig. 2(a) and 
Ref, 3]; the deep-dopant description ascribes passivation 
to the formation of a silanic bond in the vicinity of an 
^-hybridized acceptor [Fig. 2(b) and Ref. 5]. 

Much more complex is the situation at high concentra- 
tion, where unsaturated or ionized hydrogen con- 
figurations (AB, BC or HJ) are expected to rearrange 
and form saturated configurations (e.g., molecular hydro- 
gen, Hj), this process being counterbalanced by static fac- 
tors that make the incorporation of the molecules in the 
lattice difficult. 

The existence of H 2 in crystalline silicon is highly con- 
troversial: Quantum calculations show that H 2 ?a an in- 
terstitial T d site oriented in a < 100) direction is one of 



the most stable compounds in silicon and is therefore ex- 
pected to be one preferred configuration of high concen- 
tration in crystalline, nondefective silicon^ 7 however, no 
direct spectroscopic evidence for H 2 has yet been present- 
ed. On the experimental side, H 2 has been advocated in 
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FIG. 1. Suggested hydrogen configuration in silicon: (a) 
bond-centered, (b) antibonding, (c) H* , and (d) silanic. 
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FIG. 2. Suggested mechanisms of acceptor passivation in (a) 
the standard theory of shallow dopants (by a positively charged, 
bond-centered hydrogen), and (b) the deep*dopant description of 
an acceptor (by a silanic bond with an unsaturated silicon 



many cases: Hall has explained experimental diffusion 
profiles in terms of coexisting H and H 2 in local thermo- 
dynamic equilibrium; 8 Cerofolini and Ottaviani, in order 
to explain their channeling Rutherford-backscattering- 
spectrometry CRBS) data, have assumed that atomic hy- 
drogen, remaining in the crystal after low-temperature 
(77 K) ion implantation, reverts to R 2 after annealing at 
moderate temperature (400-450 'Of Stutzmann et al n 
basing their arguments on infrared and effusion data, 
reach the conclusion that Zi 2 is formed in ynneafcd p-type 
silicon; 10 but Boyce remarks that H 2 has been detected 
only in micro crystalline silicon* where it can be embed- 
ded in microcavities among micro crystals. 11 

Together with these pointlike defects, extended defects 
involving hydrogen have also been reported: namely 
(100) and (111) platelets"" 15 and bubbles. 16 

To deepen our knowledge of the state of high- 
concentration hydrogen in silicon, we have hence under- 
taken an extended experimental characterization of 
single-crystal silicon implanted at high fluence (1.6X 10 16 
cm 2 ) and low energy (15.S kcV) at room or liquid- 
nitrogen temperature with the goal of verifying under 
which annealing conditions the implanted hydrogen is 
able to produce a displacement field in the crystal. The 
techniques used in the present characterization are 
secondazy-ion-mass spectroscopy (SIMS), channeling 
RBS, elastic-recoil detection analysis (ERDAJ, and con- 
ventional and high-resolution transmission electron mi- 
croscopy (TEM and HREM). 

IL EXPERIMENT 

The complete experimental setup of our collaboration 
is described in Rc£ 17. The major difference with respect 
to that work is that all the implantations have been car- 
ried out taking care not to produce any appreciable tar- 
get heating. 

A. Sample preparation 

Samples were prepared by implanting Hj into single- 
crystal silicon. The silicon slices were 6 in. in diameter, 



Cxochralsld grown, (100) oriented, and p type, with a 
resistivity in the interval 30-50 ft cm (the corresponding 
dopant concentration is much lower than hydrogen con- 
centration, so the sample can be considered intrinsic). 
The implantations were carried out by tilting the slices by 
r with respect to the beam. 

Samples were implanted with at an energy of 31 
keV. Because of the K} fragmentation at the surface and 
low density of the coUisionaf cascades, this implantation 
is equivalent tb the implantation of atomic hydrogen at 
double fluence and an energy of 15.5 keV. Data will 
henceforth be presented in terms of atomic energy E and 
fluence*. 

The implantations were carried out either at approxi- 
mately 300 K (in an Eaton Nova 10/160 ion implanter, 
operated at a beam-current density for which there is no 
appreciable sample heating), or at approximately 77 K (in 
a Varian ion implanter operated at very low current den- 
sity, » 1 ji A/cm 2 , with an end station modified to allow 
the target to be cooled by a good thermal contact with a 
liquid-nitrogen heat sink). 

Annealing experiments were carried out by heating the 
samples in a vacuum in the temperature interval 
1Q0-8G0 # G 

• B. Measurements 

Secondary-ion mass spectrometry. Depth profiles of 
sputtered negative tons with masses of 1 amu were deter- 
mined by using a CAMECA IMS-4f ion microprobe, with 
Cs + as the primary ion at an energy of 14.5 keV, and a 
sample-chamber vacuum of about 5X10 -7 Fa. The pri- 
mary current density was adjusted to obtain a high 
sputtering rate (=3 nm/s) stable during the measurement 
to within 0.5%. These conditions allow us to reduce typ- 
ical problems of hydrogen depth profiling, mainly related 
to the contamination of the sample chamber. Knock-on 
artifacts affect the hydrogen profiles particularly in the 
distributions wines, where the hydrogen concentration is 
lower than 2X 10 Ta cm' 3 . This effect is responsible for an 
uncertainty in the depth scale smaller than the depth 
resolution (£10 nm) and for a background level increase. 
To avoid the characterization of artifacts and of diffused 
profiles, we limit our considerations around the maxima 
of the distributions. 

Elastic-recoil detection analysis. Forward scatterings 
were performed by using a 2.3-MeV He + beam (with a 
section of 0, 1 X 3 mm 2 ) impinging the sample at a glanc- 
ing angle of 75* from the sample normal The detector, 
placed 150" from the beam direction and 15 cm from the 
sample, was covered with a collimator 2 mm wide and 5 
mm long. The forward-scattered helium particles were 
stopped by a 10-pm-thick Mylar film. These data were 
used to follow the hydrogen evolution after annealing in 
the concentration region just above 10 20 cm" 3 . Com- 
parative measurements produced coincident SIMS and 
ERDA profiles within experimental uncertainties. 

Channeling Rutherjbrd-backscattering spectrometry. 
Channeling RBS analyses were performed by using a 2- 
MeV facility; the RBS spectra were obtained by imping- 
ing a 1-MeV He* beam into implanted samples, with a 
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TABLE I. Average diameter iSi and total number Gitem) of extended defects, number of atoms ar- 
ranged in extended defects Wtem), and number of displaced silicon atoms as determined from channel- 
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The sflxcon density in defects was assumed to be 10 nm" 1 . 
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scattering angle of 130* and an integrated charge in the 
range 30-100 /iC Spectra were compared with those ob- 
tained from an undamaged reference sample both in 
channeling and in random conditions. Details describing 
how the distribution of displaced atoms is extracted from 
a channeling RBS spectrum are given in the Appendix. 

Electron microscopy. Transrnission-electron-micros- 
copy observations concerning density, distribution, and 
location of extended defects were carried out with a Phi- 
lips EM 400 ST electron microscope; high-resolution in- 
vestigations of their structures were carried out with a 
JEOL 2010 electron microscope (point raoludon 0.23 
nm). 

m, EXPERIMENTAL RESULTS 

The results of interest are summarized in Figs. 3-8 and 
in Table I. 
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FIG. 3. Channeling RBS spectra of the samples implanted at 
300 K and annealed for 2 h in the temperature range 
200-400'C 



A. Formation of a displacement field 

Annealing in the range 200-400 *C induces the forma- 
tion of a strong displacement field (detected by RBS), 
quite irrespective of the target temperature during the 
implantation. Figure 3 shows the RBS spectra after 
isochronal (2 h) annealing of samples implanted at 300 K; 
no appreciable difference (apart from the as-implanted 
spectrum 9 ) was detected in the samples implanted at 77 
K. Figure 4 shows the kinetics of the formation of a dis- 
placement field in isothermal conditions (r«200"C), 
showing that this phenomenon is not saturated after 2 h 
of annealing. (In Ref. 18 it was tacitly assumed that the 
displacement field increases linearly with time; actually, 
Fig. 4 shows a complex behavior that makes the con- 
clusions reached therein concerning the activation energy 
of the process questionable.) 

Rutherford-backscattering dechanneling measurements 
(carried out at different beam energies) in samples an- 
nealed in the interval 450-500*C for 2 h show that the 
scattering centers are uncorrelated and their total 
amount per unit area is approximately 10 u cm" 2 , ir- 
respective of target temperature (77 or 300 K) during im- 
plantation. 
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FIG. 4. Kinetics of displacement-field formation (integral of 
the displacement field vs annealing time). The line has been 
drawn only to guide the eyes. 
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B. Hydrogen effiuloa and destruction 
of the displacement field 

While heat treatments at T <400°C do not produce ap- 
preciable hydrogen out-diffusion, heat treatments at TZ 
400 *C produce a fast hydrogen loss, terminating at a sat- 
uration value. This conclusion follows from the kinetics 
of hydrogen loss determined by ERD A (Fig. 5). As hy- 
drogen is gradually lost, the displacement field is progres- 
sively reduced as shown by RBS (Fig. 6). Once the dis- 
placement field has reached its maximum value (around 
400*0, heat treatments at higher temperature simultane- 
ously produce a loss of hydrogen and a reduction of the 
displacement field. The displacement field decreases with 
T approximately like the amount of hydrogen, thus sug- 
gesting that they are strictly correlated. 

The kinetics of hydrogen loss are associated with two 
regimes: the initial regime is characterized by a fast hy- 
drogen loss; the steady-state regime is characterized by a 
hydrogen concentration almost independent of the an- 
nealing duration. Interpreting this quasi-steady-state as 
an equilibrium, an effusion pseudoenthalpy LH^ can be 
determined by an Arrhenius plot of experimental data; 
this plot gives AJ5Teff-0.25-0.30 eV. 

C Extended defects 

Electron microscopy shov/s the presence of a heavily 
damaged region with depth between approximately 160 
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FIO. 6. Channeling RBS spectra of the samples implanted at 
300 K and annealed for 2 h in the temperature range 
400-800*C 
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FIO. 5. Hydrogen effusion kinetics (amount of hydrogen, 
Qu» vs annealing time f) in the temperature range 400-550 *C 
The experimental data have been fitted by the equation 
2H<rt»GH(0)expC-f/T)+GH^ where r is the lifetime of the 
effusion kinetics and is the steady-state value. 



and 260 nm, where RBS detects the presence of a strong 
displacement field both in as-implanted samples and in 
samples annealed at 400*C 

The extended defects detected in the implanted sam- 
ples are (100) and (111) platelets, often observed in 
high-flucnce hydrogen implantation, 13 and {113} plate- 
lets, usually detected in silicon-implanted layers. 19 Fig- 
ure 7 shows HREM pictures of (a) a { 100] platelet, (b) a 
{111} platelet, and (c) a (113} platelet. Through-focal 
observations in the diffraction-contrast mode have shown 
a behavior characteristic of zones with a density lower 
than that of the host matrix, therefore suggesting a possi- 
ble hydrogen content 

The shape and size of the hydrogen-related platelets 
vary with annealing temperature, ranging from a diame- 
ter of about 7 nm in as-implanted samples (300 K) to an 
average diameter of approximately 12 nm after annealing 
at 400 *C Annealing at higher temperature produces 
different defects: extrinsic defects with size increasing 
from 30 to 40 nm as T increases from 600 to 800 'C, and 
three-dimensional defects, with lower electro-optical den- 
sity, whose evolution in shape and size is understood by 
comparing the pictures shown in Fig. 8. 

Also, the in-depth distribution of extended defects 
varies with annealing temperature: it is almost uniform 
in the region 160 <x <260 nm in the samples heated up 
to 600 "C, and extends approximately over the same re- 
gion, but is bell shaped with a maximum at 190 nm in the 
samples heated at 800 'C Table I reports the most im- 
portant quantitative features of the defects observed in 
the considered samples. 



BEST AVAILABLE COPY 



4* 



HYDROGEN-RELATED COMPLEXES AS THE STRESSING . . . 



7065 



(a) 



SSI •■ > 



(b) 



Taking into account that (1) the energy behavior of 
ri c c h a n ne li ng suggests that the scattering centers are not 
correlated, €2) the amount of silicon atoms in extended 
defects and the distortion field they produce increase 
monotonically with T up to 600 # Q while the number of 
displaced atoms has a maximum of 400 and (3) for all 
T (7 £600*0 the in-depth distribution of atoms in ex- 
tended defects is almost uniform in the range 
160-260 nm, while the distributions of displaced atoms 
are quite sharp and centered on somewhat different 
depths at different temperatures, it may reasonably be ar- 
gued that the extended defects detected by electron mi- 
croscopy play a small effect on the displacement field 
detected by RBS. 
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FIG. 7. High-resolution electron micrographs of (a) a ( 100) 
platelet, (b) a {111} platelet, »"d (c)a [113] platelet 
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FIG. 8. (a) Near-<100>^rcjeenon plan view of i defective 
region in a sample annealed at 600 "C for 2 a. Through focal ob- 
servations in the dlnractton-contrast mode have shown that the 
particles axe less dense than the matrix, (b) HREM image of 
typical defects observed in samples annealed at 600 'C (c) 
Near- ( 110 > -projection cross-sectional view of defects observed 
in samples annealed at 800*0. Comparison of (b) and (c) sug- 
gests that the defect-growth mechanism is Ostwald ripening. 



IV. SOURCES OF THE DISPLACEMENT FIELD 

Since thermal treatments in the range 200-400 9 C in- 
duce the formation of a strong displacement field (Fig. 3) 
while the hydrogen profile remains substantially un- 
changed! we are forced to ascribe the displacement field 
to some hydrogen-related complexes whose nature 
changes with the annealing conditions rather than to hy- 
drogen itself. Two kinds of hydrogen-related complexes 
are identified as sources of displacement fields: one (re- 
ferred to as H 9 ) responsible for a weak displacement field, 
formed in the defect-free wing of the distribution (cen- 
tered on a depth of about 230 nm), and disappearing after 
annealing at r>200 # C; and the other (referred to as 
H„) responsible for a strong displacement field, formed 
in the defective region of maximum hydrogen concentra- 
tion after heat treatments at T> 200*0, and stable up to 
at least 6*00 

The existence of two hydrogen-related stressing species 
with different thermal behaviors was identified in previ- 
ous works of this collaboration. 30 However, due to in- 
herent difficulties in the quantitative analysis of the RBS 
spectrum, the RBS data were treated in a semiquantita- 
tive manner. Consequently, our prdiminary analyses 
may have led to erroneous conclusions. 

In this work we deepen our study by means of the new 
method for the analysis of the channeling RBS spectrum 
described in the Appendix. The analysis of the RBS spec- 
tra shown in Figs. 3 and 6 is the point of origin of the dis- 
tributions of displaced atoms reported in Figs. 9 and 10, 
respectively. The interpretation of these distributions in 
terms of a hydrogen-induced displacement field is not 
straightforward and requires that the total displaced 
atoms are subtracted by the sum of the displaced atoms 
associated with extended defects and the self-interstitials 
produced by the implantation itself. 

First of all, we recall that most of the phenomenology 
observed by RBS is not due to extended defects. Indeed, 
the correlation of the scattering centers (uncorrelated 
versus correlated), the annealing behavior (the RBS signal 
is at a maximum at 400*C, while extended defects in- 
crease monotonically up to 600 *Q, and the form of the 
distribution (centered on a depth between 200 and 250 
nm according to T with a full width at half-height of 50 
nm versus a uniform distribution in the range 160-260 
nm] give evidence that extended defects bring a small (but 
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not negligible) contribution to the displacement field 
detected by RfiS; However, since in the sample annealed 
at 800 'C hydrogen is no longer present and the in-depth 
distribution of the displaced atoms detected by channel- 
ing RBS is the same as the in-depth distribution of ex- 
tended defects detected by TEM, it may reasonably be ar- 
gued that in this situation all the displaced atoms are due 
to extended defects. This argument leads to the con- 
clusion that each atom in the extended defect is responsi- 
ble for approximately 20 equivalent RBS-active atoms 
(see Table D. In all other situations in which hydrogen is 
present the in-depth distribution of displaced atoms pro- 
duced by extended defects is assumed to he the same (flat 
up to 600 *C and bell shaped at 800 *C) as determined by 
TEM analysis multiplied by the above motor 20. 

The self-interstitial distribution can be calculated by 
existing simulators On our case, TRIM, denoting "trans- 
port of ions in matter** 21 ) by adjusting the electronic stop- 
ping (by approximately \0%) to match the experimental 
hydrogen distribution and the threshold displacement en- 
ergy to match the displaced -a torn distribution in the re- 
gion where the effect of the implanted impurity is negligi- 
ble (in our case, in the depth region x 5 120 nm where 
neither H* nor H„ are present). 

The analysis of the RBS spectra of the samples im- 
planted at 300 Kin the depth region x % 120 nm suggests 
that the as-implanted self-interstitial distribution is un- 
changed after annealing at r£300*C and is appreciably 
annealed out for 7'fc400 , C This conclusion is consistent 
with the thermodynamic analysis of Roorda et al. show- 
ing that the recovery of crystalline silicon damaged by 
the implantation of high -energy, fight ions takes place 
quite uniformly in the temperature range 100-400 'C, the 



released energy being approximately I eV per displaced 
atom. 22 This' value results from the weighted average of 
the energy excess stored in Frenkel pairs [tentatively, in 
the range 4*8 eV/pair (Rcf. 23)}, amorphous islands 
formed by high-energy silicon recoils 24 (approximately 
0.13 eV/atom (Rcf. 25)], or agglomerates of point defects 
(eg., divacancies 26 ). Since the recovery kinetics of 
separated amorphous islands, formed by arsenic implan- 
tation (the . ion for which data are available with the 
closest mass to that of silicon, 27 is in the range 
100- 150 *C, while divacancies anneal out around 
200 "C, 24 the following attribution seems realistic: 

(a) 100-1 50 *C temperature range for the annealing of 
amorphous islands. 

(b) 150-250 'C, temperature range for the annealing of 
divacancies. 

(c) 250-400 'C, temperature range for the annealing of 
selfSnterstitials. 

Remembering that in our experiment' the implantation 
energy was not sufficient to form amorphous islands and 
that divacancies are, to the first order, RBS-inactive, it 
may reasonably be argued that the RBS signal due to 
self-inters titials disappears in the temperature range 
250-400 "C, in substantial agreement with our finding. It 
is also possible that in our case the system is partially 
protected against Frenkel-pair recombination by vacancy 
or interstitial reaction with the atomic hydrogen. 
• Subtracting the remaining inters titials and the dis- 
placement field due to extended defects from the-Eotal.dis- 
placed atoms, the hydrogen-induced distributions are 
eventually obtained. An example showing the various 




1.0 



SO 100 150 200 250 300 
Depth (nm) 

FIG. 9. In-depth distributions of displaced atoms of samples 
annealed at a temperature below 400 *C for 2 h; beating at 
350 "C produces the same displaced-atom distribution as heating 
at 400 X. 
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FIG. 10. In-depth distributions or displaced atoms of samples 
annealed in the temperature range 400-800 *C for 2 h. 
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FIG. 11. Separation of the various contributions (hydrogen 
induced, or due to self-interstftlals or extended defects) to the 
displacement field in the sample * "tinted at 200 "C 



contributions to the experimental displacement field in 
the sample annealed «t 200 "C is shown in Fig. 1 1. 

The above analysis confirms the original qualitative 
conclusion relating the displacement field to hydrogen- 
related complexes: A displacement field is already 
present in as-implanted conditions (300 K) in the bulk tail 
of the hydrogen distribution (x £250 nm) and remains 
unchanged after annealing at 150 *C; heating at a temper- 
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FIO. 12. Graphs showing the evolution vs annealing temper- 
ature of the amount of hydrogen-displaced silicon atoms and or 
their most probable depth. The lines have been drawn only to 
guide the eyes. 



ature higher than 200 *C (but not greater than 400*0 
produces an increase in intensity and a gradual shift to- 
ward the surface of the displacement field; heating at 
higher temperature produces a gradual loss of hydrogen 
without further change of shape or depth distribution 
(Fig. 12). 



V, STRUCTURE OP THE HYDROGEN-RELATED 
COMPLEXES 

The analysis described in the preceding sections has led 
to the identification of two stressing species, referred to 
as H. and H.„ irrespective of their actual chemical na- 
ture. 

In view of its depth location and annealing behavior, 
the identification of H # as molecular hydrogen, H 2 , seems 
plausible; of course, this attribution has mainly a conven- 
tional character, and in the absence of direct spectroscop- 
ic evidence it must be considered a putative stipulation. 

In this section we speculate on the nature of the H,. 
center. The driving force leading to the formation of the 
H 2 molecule starting from silanic, BC or AB 
configurations is the electronic energy gained in its syn- 
thesis (see, for instance, the estimate of Ref. 20), which 
greatly counterbalances the modest energy required to 
sustain the weak displacement field produced by the mol- 
ecule. 

What is obscure in this picture is which driving force is 
able to transform H 2 molecules, each producing a weak 
displacement field, in the complex H„, which .requires 
an appreciable energy to sustain the strong displacement 
field it originates. Molecular hydrogen in interstitial sites 
is, however, in a strained configuration, as confirmed by 
the calculation of its bond length (0;086 nm, 15% larger 
than in vacuum) and energy (approximately 2 eV, versus 
4.75 eV of the free molecule).* 1 An H 2 molecule in a va- 
cancy or in larger voids is therefore expected to be more 
stable, and its formation is allowed because the energy re- 
quired to dissociate two silanic bonds is counterbalanced 
by the molecule formation and the saturation of the re- 
sulting dangling bonds. 20 Since the H„ complex is 
formed in the zone of maximum hydrogen concentration, 
where intrinsic defects are presumably fully decorated, it 
is therefore not unrealistic to assume that H,, is formed 
by several H 2 molecules in defects of intrinsic nature, 
such as vacancies, divacandes, intrinsic stacking faults, 
etc. 

The location of the H„ complex in the zone of max- 
imum radiation damage, together with the fact that after 
annealing at 7*>400*C the displacement field decreases 
approximately in proportion to the total amount of hy- 
drogen, might suggest that the H„ complex is stabilized 
by its interaction with silicon extended defects formed 
during annealing. The decorated defects so formed could 
eventually evolve to form the (100) and (111) platelets 
observed in Rets. 12 and 13 by transmission electron mi- 
croscopy. Though this model is partially sustained by 
TEM evidence that- (100 J, (111), and {113} defects ob- 
tained after annealing at 400 *C are less dense than bulk 
silicon (thus suggesting their possible decoration with hy- 
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drogen), it does not explain the RBS behavior for 

r<4oo*a 

More plausible is the following mechanism, for which 
H„ is seen as a set of separated molecules fixed in prefer- 
ential sites (vacancy clusters) by steric factors. 

(1) In as-implanted conditions at 300 K v atomic hydrogen 
(in BC or interstitial positions) diffuses until it reacts with 
another hydrogen atom to form a hydrogen molecule, 

H+H— H 2 , 

which rests in the position where it was formed. This 
process is expected, to take place in the region where va- 
cancies to) either are not formed or are in insufficient con- 
centration to bind hydrogen with silanic bonds and is 
completed at T£ 200 *C; this argument suggests that H 2 
is formed starting from x Z2QO run to a depth limited by 
atomic hydrogen diffusivity. 

(2) Annealing at Tj£200*C produces decomposition of 
silanic bonds 21 of fully decorated vacancies (denoted u- 
H 4 ) and the formation of the H*« complex, 

in the region where close, fully decorated vacancies are 
present (150-200 nm). The hydrogen-related defect so ' 
produced is formed by n/2 vacancies embedding n H 2 
molecules (¥L m ^o 9/1 • • • nHj) and is thought of as im- 
mobile because the close vacancies form a preferential 
site where hydrogen can stay. 

This mechanism (which requires the existence of intrinsic 
defects, each bonding several hydrogen atoms, evidence 
of which has recently been given 29 ' 38 ) also explains why 
H,. is formed in the zone with simultaneously high hy- 
drogen and vacancy concentrations. 

The size, of the hydrogen-related complex H„ can be 
estimated from RBS dechanneling data by assuming that 
after prolonged annealed at 7*>400*C practically all hy- 
drogen is embedded in separated complexes and that each 
of them displaces approximately 10 3 atoms from their lat- 
tice positions (this estimate is taken from the thermo- 
dynamic analysis of Kef. 23, Chap. 7); this estimate gives 
that the H„ complex is formed by approximately four 
H 3 molecules (this gives a better specification of the nano- 
cavity at 400*C, H„, «d 2 • • • 4H 2 ). 

Assuming that each H„ complex is actually formed 
by n molecules in a supercompressed gas phase at a densi- 
ty of the order of (0.5-1 )X 10 23 cm (corresponding to 
one to two H 2 molecules per vacancy), the pressure exert- 
ed by the gas on the silicon is of the order of 10 B Pa and 
the corresponding lattice strain, =*10~ 3 , is of the same 
order as that observed by this collaboration by multiple- 
crystal x-ray diffraction. 

The nanocavities in which the complexes are em- 
bedded become closer and closer as fluence increases and 
eventually may coalesce to form voids containing hydro- 
gen bubbles, as reported in the literature. 15 - 16 We have 
also been convinced of the formation of cavities by the 
observation of octahedral defects (which are the equilibri- 
um shapes of cavities in silicon) after annealing at 800 *C 



[see Fig. 8(c)]. 

This model also explains why the system prefers to ar- 
range in a configuration of higher stress — indeed, the 
driving force is the energy gained in relaxing the H-H 
bond from a length of 0.086 nm (bond energy of =2 eV) 
in its interstitial position to a length of 0.074 nm (bond 
energy of 4.75 eV) in vacancies or larger cavities. 

^ VX CONCLUSIONS 

A combined experimental analysis of the Si Ji system 
obtained by high-fluence ion implantation in (100) single- 
crystal silicon at room or liquid-hydrogen temperature, 
based on channeling RBS, conventional and high- 
resolution TEM, ERDA, and SIMS profiling, has allowed 
us to identify the following important items. 

(1) Two hydrogen-related complexes, H, and H„ are 
the major sources of the displacement field in the crystal. 

(2) H, is observed in the defect-free wing of the hydro- 
gen distribution in as-implanted samples and is destroyed 
by annealing at r>200*Q this hydrogen-related center 
is attributed to isolated H 2 molecules. 

(3) H„ is formed after annealing at 7>200*C in the 
silicon region where radiation damage and hydrogen con- 
centration are both high; this hydrogen-related center is 
stable up to at least 600*C and is attributed to small bub- 
bles of H 2 molecules. 
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APPENDIX: ANALYSIS OP RBS SPECTRA 

In this appendix we shall assume that the energy axis 
of RBS spectra has already been transformed in a depth 
axis and that the backscattering yield Y as a function of 
depth x is experimentally known. Together with the 
function Y{x) for the considered sample, two other RBS 
profiles are assumed known: Y ttzdmtA {x) 9 the RBS yield 
for an undamaged crystal, and Y na * m {x) 9 for the ran- 
domly oriented (or amorphous) silicon. 

Channeling RBS spectra are often used in a qualitative 
fashion, since their quantitative interpretation requires 
the exact specification of the lattice location of the back- 
scattering centers 17 and a knowledge of how multiple 
scatterings, in turn responsible for dechanneling, affect 
the baseline. This last point is especially important in sit- 
uations (like those considered in this work) in which two 
peaks of comparable intensities are suspected to be 
present and separated by a distance (=40 nm) just slight- 
ly higher than the depth resolution (=* 20-30 rim) of the 
RBS technique. The peaks in Figs. 3 and 6 indeed seem 
centered on different depths in samples annealed at 
different temperatures, but this conclusion might only be 
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apparent because the (unknown) baseline B{x) affects the 
estimate of the direct backscattering yield D(x) 
[Z>(x)=r(x)-£(x)]. This effect is especially important 
in all cases dike ours) in which B (x) is of the same order 
asi)(x). 

We have hence undertaken the development of an 
objective method for constructing the baseline on serious 
physical bases; this method seems to us quite general. 

The basic assumption upon which our method is con- 
structed is the following: 

(A) The concentration of backscattering centers is pro- 
portional to the concentration of dechanneling centers. 

This hypothesis allows the baseline B {x) to be construct- 
ed: 

With reference to the ideal situation of Fig. 13 (which 
is representative of a wide class of actual cases), first of all 
identify the depth region (x Q ,x P ) in which scattering 
centers are present; in principle, Xq can coincide with the 
surface. In the region (x Qt x P ) the baseline increases 
gradually from BW-Y^^ix) for x<x Q to 
B{x)=YU) for x P <x. The analysis cannot proceed fur- 
ther without approximating the baseline in this region by 
means of an analytical function; a linear approximation, 
£(x)- y>+m(x— x P ) [where m is the tangent to the ex- 
perimental curve Y(x) in a right neighbor of P], is partic- 
ularly useful. 

The gradual increase of the baseline in the interval 
(xq,x p ), physically due to the loss of collimation of the 
impinging beam as a greater and greater number of 
scattering centers are met, can be described by an equa- 
tion such as 

B(x)~[l-Mx)]Y' l ^x)+k(x)[Yp+m(x--Xp)] . 

(Ai) 

The dependence of k on depth x is specified by the as- 
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FIG. 13. Procedure for the calculation of the baseline in the 
RBS spectrum On principle, x Q can coincide with the surface, 
and the baseline is manifestly stable with respect to the choice 
of x, provided x f £ 300 nm). 



sumed model of dechanneling. Assuming that there is no 
additional dechanneling (compared to the undamaged 
sample) until all of the damaged layer has been crossed by 
the beam, one has 



Vx: x 0 <x<xp— &(x)— 0 



CA2) 



(this choice produces a lower limit to the baseline); taking 
Hx) such that 



*(x)« 



r^x) forO<x<x ff , 

Y Q +( Y P - Y Q )(x -x ff )/(x P -x Q ) for x Q <x <x P , 
Y(x) otherwise , 



(A3) 



one encounters the frequently used "linear dechanneling 
approximation." 
We quantitatively specify assumption (A) by assuming 

X(x)« f*D(x)dx • 
*Q 

The constant of proportionality is calculated by speci- 
alizing this relationship to the endpoint P t where 
X(x P )»l,sothat 

Ux)=rDU)dx/r F D<x)dx = ArD{x)dx , 

(A4) 

with 



Inserting (A4) into (Al), one gets the following integral 
equation: 

B(x)^Y^ix) [l-Af^[Yix)-B{x)]dx} 

+ [Y P +m{x-x P )]A [ x lYlx)-B(x)]dx 

which can be solved iteratively starting, for instance, 
from the lower bound £ 0 (x)= Y^^x). 

Once the baseline is calculated, the profile of direct 
backscattering yield is immediately calculated. This 
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profile is eventually transformed into a concentration 
profile by "normalization to random," Lc, by assuming 
that the concentration of backseat tering centers is given 
by 

Cix)^n Si D{x)/D nttaom {x) , (A5) 

where D aadem {x) is the difference of the random signal 
(coming from a target where all atoms are displaced) 
minus the baseline (this assumption is realistic for* base- 
lines appreciably lower than the random signal). 

The method hitherto described indeed seems objective, 
as it does not depend on the points Q and P (while point 
Q is 'unequivocally determined from the experimental 



data, point P is, in principle, arbitrary). While both 
choices (A2) and (A3) have solutions that depend 
markedly on the somewhat arbitrary choice of P> the 
method proposed here is very stable with respect to this 
choice in a wide depth interval, which makes us confident 
of the robust physical basis of our computed displace- 
ment field. 

It must be noted, however, that in Eq. IA5) a factor 
0(1) is ignored because the backscattering yield of a dis- 
placed atom (however large is its displacement) has been 
assumed to be the same as in a randomly oriented crystal 
Moreover, if two different displacement fields are present 
(as, for instance, in the case of hydrogen), two factors, in 
principle different and both O ( 1 ), are ignored. 
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